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Abstract 

We study the discovery potential of future beta decay experiments on searches 
for the neutrino mass in the sub-eV range, and, in particular, KATRIN experiment 
with sensitivity m > 0.3 eV. Effects of neutrino mass and mixing on the beta decay 
spectrum in the neutrino schemes which explain the solar and atmospheric neutrino 
data are discussed. The schemes which lead to observable effects contain one or two 
sets of quasi-degenerate states. Future beta decay measurements will allow to check 
the three neutrino scheme with mass degeneracy, moreover, the possibility appears to 
measure the CP-violating Majorana phase. Effects in the four neutrino schemes which 
can also explain the LSND data are strongly restricted by the results of Bugey and 
CHOOZ oscillation experiments: Apart from bending of the spectrum and the shift 
of the end point one expects appearance of small kink of (< 2%) size or suppressed 
tail after bending of the spectrum with rate below 2 % of the expected rate for zero 
neutrino mass. We consider possible implications of future beta decay experiments for 
the neutrino mass spectrum, the determination of the absolute scale of neutrino mass 
and for establishing the nature of neutrinos. We show that beta decay measurements 
in combination with data from the oscillation and double beta decay experiments will 
allow to establish the structure of the scheme (hierarchical or non-hierarchical), the 
type of the hierarchy or ordering of states (normal or inverted) and to measure the 
relative CP-violating phase in the solar pair of states. 
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1 Introduction 



The reconstruction of the neutrino mass spectrum is one of the fundamental problems of 
particle physics. The program includes the determination of the number of mass eigenstates, 
and of the values of masses, mixing parameters and CP-violating phases. 

At present, the evidence for non-zero neutrino mass follows from oscillation experiments 
which allow to measure the mixing parameters \U a j\, the mass squared differences and, 
in principle, the so called Dirac CP-violating phases. However, the absolute values of the 
neutrino masses cannot be determined. From the oscillation experiments one can only extract 
a lower bound on the absolute value of neutrino mass. Obviously, for a given Am 2 , at least 
one of the mass eigenvalues should satisfy inequality: 

rrii > \J\Krn?\. 

Thus, the oscillation interpretation of the atmospheric neutrino data |l| gives the bound: 

m 3 > y/Am 2 atm ~ (0.04 - 0.07) eV. 

Clearly, without knowledge of the absolute values of neutrino masses our picture of Nature 
at quark-lepton level will be incomplete. The knowledge of absolute values of neutrino masses 
is crucial for understanding the origin of the fermion masses in general, the quark-lepton 
symmetry and unification. The determination of the absolute mass scale of neutrinos is 
at least as important as the determination of other fundamental parameters such as the 
CP-violating phases and the mixing angles. Actually, it may have even more significant 
and straightforward implications for the fundamental theory. It is the absolute mass which 
determines the scale of new physics. 

The absolute values of masses have crucial implications for astrophysics and cosmology, 
in particular, for structure formation in the Universe. In fact, the recent analysis of the latest 
CMB data (including BOOMERanG, DASI, Maxima and CBI), both alone and jointly with 
other cosmological data (e.g., galaxy clustering and the Lyman Alpha Forest) shows that [[| 

m u < 2.2 eV, (1) 

for a single neutrino in eV range. Future observations can improve this bound. The Planck 
experiment will be sensitive to neutrino masses down to m v ~ 1 eV [[|. However, the 
cosmological data may not be conclusive. Even if some effects are found, it will be difficult 
to identify their origin. Modification of the original spectrum of the density fluctuations 
can mimic to some extent the neutrino mass effect. If no distortion is observed in the 
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spectrum, one can put an upper bound on the neutrino mass assuming, however, that there 
is no conspiracy which leads to cancellation of different effects Therefore independent 
measurements of the neutrino mass are needed and their results will be used in the analysis 
of the cosmological data as an input deduced from particle physics. 

Several methods have been proposed to determine neutrino masses by using the supernova 
neutrino data. One method is based on searches for the energy ordering of events which 
has, however, rather low sensitivity ||. The limits on the mass can be also obtained from 
observations of sharp time structures in the signals. It was suggested to study the time 
distribution of detected neutrino events emitted from supernova which entails to black hole 
formation ||. By this method Super-Kamiokande can measure values of the v e mass down 
to 1.8 eV and SNO can put an upper bound 20 eV on the and v T masses . (Clearly this 
bound on the and v T masses is much weaker than bounds implied by combined analysis 
of the solar and atmospheric neutrino data and direct measurements of the u e mass.) In 
this case one can check the still non-excluded possibility in which the solar neutrino problem 
is solved by the oscillations to sterile neutrino and the masses of z/ M and v T are in 20 eV 
range. (Such neutrinos should be unstable in cosmological time.) The absolute values of the 
neutrino masses can be determined in the assumption that the cosmic rays with energies 
above the GZK cutoff are produced in annihilation of the ultra-high energy neutrinos with 
the cosmological relic neutrinos |7|, ||. The analysis of the observed energy spectrum of 
cosmic rays above 10 20 eV gives the mass m u = (1.5 — 3.6) eV, if the power-like part of the 
ultra-high energy cosmic rays spectrum is produced in Galactic halo, and m u = (0.12 — 0.46) 
eV, if this part has the extragalactic origin [|10] . 



Neutrinoless double beta decay (2/30z/) searches are sensitive to the Majorana mass of the 
electron neutrino. However, in the presence of mixing the situation can be rather compli- 
cated: The effective Majorana mass of v e relevant for the 2/50z/-decay, combination 
of mass eigenvalues and mixing parameters given by 



(2) 



From this expression it is easy to find that if the 2/30^-decay is discovered with the rate 
which corresponds to m ee , at least one of the mass eigenvalues should satisfy the inequality 
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rrii > — — , (3) 
n 

where n is the number of neutrino mass eigenstates that mix in the electron neutrino. This 
bound is based on the assumption that exchange of the light Majorana neutrinos is the only 
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mechanism of the 2/30z/-decay and all other possible contributions are absent or negligible. 
Another uncertainty is related to n. We know only the lower bound: n > 3. 

The best present bound on the 2/30z/-decay obtained by Heidelberg-Moscow group gives 



m ee < 0.34 (0.26) eV, 90 % (68%) C.L.. (4) 

This bound, however, does not include systematic errors related to nuclear matrix elements 

A series of new experiments is planned with increasing sensitivity to m ee : CUORICINO 
0, CUORE (m ee ~ 0.1 eV) (TJ, MOON (m ee ~ 0.03 eV) fTJ and GENIUS (m ee ~ 0.002 



eV) 0. 

Although the knowledge of m ee provides information on the mass spectrum independent 
of Am 2 's, from m ee one cannot infer the absolute values of neutrino masses without additional 
assumptions. Since in general the mixing elements are complex there may be a strong 
cancellation in the sum (|2|). Moreover, to induce the 2/30U decay, v e must be a Majorana 
particle. 

The information about the absolute values of the masses can be extracted from kinematic 
studies of reactions in which a neutrino or an anti- neutrino is involved (e.g., beta decays or 
lepton capture). The most sensitive method for this purpose is the study of the electron 
spectrum in the tritium decay: 

3 H ^ 3 He + e"+ u e . (5) 
In absence of mixing, the energy spectrum of e~ in (||) is described by 

^ = R(E)[(E - E f - ml]*, (6) 



see, e.g., flT|| ) where E is the energy of electron, E is the total decay energy and R(E) is 



R(E) = G 2 FW ^ cos 2 9 C \M\ 2 F{Z, E)pE(E - E) . (7) 



a m^-independent function given by 

4 

2tt 3 

Here Gf is the Fermi constant, p is the momentum of the electron, 9c is the Cabibbo angle 
and M is the nuclear matrix element. F(Z, E) is a smooth function of energy which describes 
the interaction of the produced electron in final state. Both M and F(Z, E) are independent 
§In what follows we will use the bound (|J) in our estimations for definiteness. At the same time, one 

should keep in mind that due to uncertanties of nuclear matrix element the values of m ee up to ~ 0.5 eV 

can not be excluded. 



4 



of m„, and the dependence of the spectrum shape on m v follows from the phase volume factor 
only. The bound on neutrino mass imposed by the shape of the spectrum is independent of 
whether neutrino is a Majorana or a Dirac particle. 

The best present bound on the electron neutrino mass, (obtained in the assumption of 
no mixing) is given by Mainz tritium beta decay experiment |L8] : 



m Ue < 2.2 eV (95% C.L.). (8) 

Analysis of the Troitsk results leads to the "conditional" (after subtraction of the excess of 
events near the end point) bound |19] 



m Ue < 2.5 eV (95% C.L.). (9) 

The present spectrometers are unable to improve the bounds (|| |^) substantially. Further 
operation of Mainz experiment may allow to reduce the limit down to 2 eV. In this connection 
a new experimental project, KATRIN, is under consideration with an estimated sensitivity 
limit HI 

m Ue ~ 0.3 eV. (10) 
In the case of negative result fromthe KATRIN searches one can get after three years of 



operation the bound m Ue < 0.35 (0.40) eV at 90 % (95 %) C.L. |0 . 

Note that with this bound KATRIN experiment can explore the range of neutrino mass 
which is relevant for the Z-burst explanation of the cosmic ray with super-GZK energies 0. 

The aim of this paper is to study the discovery potential of the next generation tritium 
beta decay experiments with sensitivity in the sub-eV range and in particular, KATRIN 
experiment. We consider the effects of neutrino mass and mixing on the /3-decay spectrum 
expected for specific neutrino schemes. We describe the three-neutrino schemes which are 
elaborated to explain the data on the solar and atmospheric neutrinos as well as the four- 
neutrino schemes which accommodate also the LSND result. We study the bounds that the 
present and forthcoming 2/30^-decay searches, as well as the oscillation experiments can put 
on possible tritium decay results. We also consider the implications of future beta decay 
measurements for the identification of the neutrino mass spectrum. 

The paper is organized as follows. In section 2 we give a general description of the effect 
of massive neutrinos on the beta decay spectrum in the presence of mixing. In section 3 
the three-neutrino schemes are explored. In section 4 we present a general discussion of 
predictions for the beta decay in the four-neutrino schemes which explain the LSND result. 
We emphasize the importance of the bounds on the beta decay parameters imposed by Bugey 
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and CHOOZ experiments. In section 5 we study the properties of the beta decay in the 
hierarchical four- neutrino schemes. In section 6, the non- hierarchical four- neutrino schemes 
are considered. In section 7 we summarize the role that future beta decay measurements will 
play in the reconstruction of the neutrino mass spectrum. Conclusions are given in section 
8. 



2 Neutrino mixing and beta-decay. The effects of de- 
generate states 

In presence of mixing, the electron neutrino is a combination of mass eigenstates Vi with 
masses mf v e — J2i U ei i>i. So that, instead of ©, the expression for the spectrum is given 
by 

^ = R(E) £ \U ei \ 2 [(E - Ef - m^e(E -E- m,), (11) 
dE i 

where R(E) is defined in (0). The step function, Q(Eq — E — m^), reflects the fact that a 
given neutrino can be produced if the available energy is larger than its mass. According to 
eq. (|TT| ) the presence of mixing leads to distortion of the spectrum which consists of ^ 

(a) the kinks at the electron energy = E ~ E$ — rrii whose sizes are determined by 
\U -I 2 - 

(b) the shift of the end point to E ep = Eq — mi, where m± is the lightest mass in the 
neutrino mass spectrum. The electron energy spectrum bends at E ~ E ep . 

So, in general the effect of mixed massive neutrinos on the spectrum cannot be described 
by just one parameter. In particular, for the three-neutrino scheme, five independent pa- 
rameters are involved: two mixing parameters and three masses. 

Substantial simplification, however, occurs in the schemes which explain the solar and 
atmospheric neutrino data and have the states with absolute values of masses in the range 
of sensitivity dlQ|) . The simplification appears due to existence of sets of quasi-degenerate 
states. Indeed, in these schemes there should be eigenstates with mass squared differences 
Amg < 2 • 10~ 4 eV 2 and Am 2 tm ~3 • 10~ 3 eV 2 . If the neutrino masses, m i} are larger than 
0.3 eV fllOD , the mass differences 

A,~^ (12) 
2m 

^In what follows we will use the terminology elaborated for the ideal Kurie plot without background. 
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turn out to be smaller than 5 x 10~ 3 eV. Moreover, Am/m ~ Am 2 /2m 2 <C 1, that is, 
the states are strongly degenerate. Since the detectors cannot resolve such a small mass 
split, different masses will entail just to one visible kink with certain effective mass and 
mixing parameter. As a consequence, the number of relevant parameters which describe the 
distortion of the beta spectrum is reduced to one or three, depending on the type of the 
scheme (see sections 3-6). 

In the Ref. |ZT| it has been shown that for energies E Ve ^> m n the distortion of the elec- 



tron energy spectrum in the /3-decay due to non-zero neutrino mass and mixing is determined 
by the effective mass 



L-ii \ u ei\ 

However, the highest sensitivity to the mass of V{ appears in the energy range close to the 



Z^i \ u ei\ 



end point where E v ~ m« and therefore the approximation used in [21] to introduce m 



does not work. In what follows we show that still it is possible to use the mass parameter 



(|13|) for a set of quasi-degenerate states. 

In general, the neutrino mass spectrum can have one or more sets of quasi-degenerate 
states. Let us consider one such a set which contains n states, Vj, j — i, i + 1, i + n — 1 
with Arriji <C rrij. We define AE as the smallest energy interval that the spectrometer can 
resolve. (Note that AE may be smaller than the width of resolution function, and the latter 
is about 1 eV in KATRIN experiment.) We assume that Amy <C AE. 

Let us introduce the coupling of this set of the states with the electron neutrino as 

Pe^El^f, (14) 

j 

where j runs over the states in the set. We will show that the observable effect of such a 
set on the beta spectrum can be described by p e and the effective mass mp which can be 
introduced in the following way. Let us consider the interval AE in the region of the highest 
sensitivity to the neutrino mass, that is, the interval of the electron energies 

{E Q - mi - AE) - {E - rm) , (15) 

where m 8 is the mass of the lightest state in the set. The number of events in this interval, 
An, is given by the integral 

An = / a —dE. (16) 

J E -m,i-AE dry 

We will define the effective mass in such a way that the number of events calculated 
for the approximate spectrum with single mixing parameter p e and mass mp, An(p e , mp), 
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reproduces, with high precision, the number of events calculated for exact neutrino mass 
and mixing spectrum An(U e j,rrij). That is, 

R = An(p e ,mp) - An(U ej ,mj) ^ 1 
~~ An{p e ,mp) 

Expanding An (see Appendix) in powers of Arrij/AE <C 1, where 

Arrij = rrij — mp , (18) 

we obtain: 

««El^l 2 ^ + 0((^)). (19) 
It is easy to see that the first term vanishes if 

mp = ^>i m i\ U °i\\ (20 ) 

Pe 

So, for this value of mp, R is of the order of (Arrij / AE) 2 . Note that if we set mp to be equal 
to the value of any mass from the set or the average mass, the difference of the number of 
events would be of the order of Am/ AE. If AE is relatively small, this correction may be 
significant. The expression (|20|) is similar to ([13]), but in (|20"D j runs over a quasi-degenerate 
set (not over all the states.) Moreover, provided that Am <C AE, the approximation works 
for all energies. 

In reality the background should be taken into account. However it is easy to see that if 
the change of the background with energy in the interval AE is negligible, our analysis will 
be valid in the presence of the background, too. 

If the scheme contains more than one set of quasi-degenerate states with the correspond- 
ing effective masses m q p and mixing parameters pi, the observable spectrum can be described 
by the following expression 

^ = R(E)J2pI [(E -Ef-(ml ff fME Q -E-ml ff ), (21) 

where q runs over the sets. Each set of quasi-degenerate states will produce a single kink at 
the electron energy E q ~ Eq — mt with the size of the kink determined by p\. The set with 
the lightest masses leads to bending of spectrum and the shift of the end point. 

3 Three neutrino scheme 

Let us consider the three-neutrino schemes which explain the solar and atmospheric neutrino 
results. In the case of mass hierarchy, m\ <C m 2 <C m 3 , the largest mass, m 3 ~ J Aml tm = 



(4 — 7) x 10~ 2 eV, is too small to produce any observable effect in the planned tritium decay 
experiments (see eq. (|10D). 

If m 3 is in the sensitivity range of KATRIN experiment (777.3 > 0.3 eV), the mass spectrum 
should be quasi-degenerate. Indeed, 



Am 3i _ Am^, 
m 3 2ml 



< 0.03. 



Moreover, from the unitarity condition we get the coupling parameter 

Pe= E \Uej\ 2 = l. (22) 
3=1,2,3 

Therefore the effect of non-zero neutrino masses and mixing on the /3-decay spectrum is 
characterized by unique parameter - the effective mass 

m P = E m j\Uej\ 2 - m 3 . (23) 

3=1,2,3 

Correspondingly, the distortion of the /?-decay spectrum consists of a bending of the spectrum 
and shift of the end point determined by mp {Eq — > Eo—mp), as in the case of v e with definite 
mass and without mixing. Let us consider the bounds on mp imposed by the 2/50i/-decay 
searches and the oscillation experiments. (The 2/?0^-decay in schemes with three degenerate 
neutrinos has been extensively discussed before [p2j| , p3| ). Assuming that neutrinos are 
Majorana particles we get from (fj) and (|20| ) the relation between the effective masses in the 
beta decay and the double beta decay: 

m ee ~ m p \\U el \ 2 + e l ^\U e2 \ 2 + e^\U e3 \% (24) 

where 02 and 03 are the relative CP-violating phases of the contributions from the second 
and the third mass eigenstates. 



According to the CHOOZ bound [24] which is confirmed by the slightly weaker bound 
obtained in Palo Verde experiment one of the squared mixing elements (let us take \U e3 \ 2 
for definiteness) must be smaller than 0.05. The other two elements are basically determined 
by the mixing angle 6 & responsible for the solution of the solar neutrino problem, so that 
the eq. (|4]) can be rewritten as 

m ee = m p \(l- \U e3 \ 2 )(cos 2 6 Q + e^ 2 sin 2 9 e ) + e^\U e3 \ 2 \ . (25) 

From this equation we find the following bounds on the beta decay mass (see also fTLj| ): 

m ee <mp< — — — — 77-^:, 26 

cos 20© (1 - \U e3 \ 2 ) - \U e3 \ 2 \ 
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where the upper bound corresponds to the maximal cancellation of the different terms in 



(H). 

The bounds ( P6f) are shown in fig. |l[ (See also [p2].) The following comments are in 
order: 

1) For zero value of U e % the weakest bound on mp from the double beta decay appears 
at maximal mixing: tan 2 # Q = 1. For non-zero U e z the points of the weakest bound shift to 
tan 2 8 & ~ 1 ± 2|[/e3| 2 . In the vicinity of these points, the upper bound on mp is given by the 
present beta decay result (see eq. (|[)). 

2) Taking the best fit values of 9 & from the various large mixing solutions of the solar 



neutrino problem [£6] we find from eq. (|26|) the following bounds: 

f 0.67 -0.74 eV LMA 
mp < < 1.6 - 2.2 eV LOW , (27) 
{ 1.0-1.3 eV VAC 

where we have used m ee <0.34 eV (^) and the two numbers in each line correspond to |C/ e 3| 2 = 
and 0.05, respectively. Note that already existing data on the 2/30^-decay give bounds (at 
the best fit points) which are stronger than the present bound from direct measurement. 

Moreover, for m ee < 0.07 eV which can be achieved already by CUORE experiment 
[14]], the bound on mp from 2/30z/-decay in the LMA preferable region of tan 2 9 is below the 
sensitivity of KATRIN experiment. Therefore, the positive result of KATRIN experiment 
(and identification of the LMA solution of the solar neutrino problem) will lead to exclusion 
of such a 3^-scheme. 

3) For the SMA solution of the solar neutrino problem we get mp ~ m ee , and conse- 
quently, according to the bound (||): mp < 0.34 eV. So, the expected range of rap only 
marginally overlaps with the KATRIN sensitivity region. Thus, if the SMA solution is iden- 
tified and the LSND result is not confirmed, favoring the three neutrino scheme, the chance 
for observation of the beta spectrum distortion in KATRIN experiment is rather small. 

4) A positive signal in the 2/50z/-decay searches will have important implications for the 
tritium decay measurements: 

a) . According to (0), it gives a lower bound on mp independently of the solution of the 
solar neutrino problem: m ee < mp. 

b) . If the values of m ee , mp and \U e3 \ 2 are measured, we will be able to determine the 



CP-violating phase (fi 2 in (|25|) : 

1 



. 2 
sin 



2 sin^ 26 r. 



mp J \ \ mp ) mp 



10 



sin 2 29 







1 _ ( m eeV 



where (±) sign of the last term reflects an uncertainty due to the phase of (03 in (]24])). 

c). If mp turns out to be smaller than m ee , we will conclude that there are some additional 
contributions to the 2/30i/-decay unrelated to the Majorana neutrino mass. 

If future (3 decay measurements with sensitivity flTD| ) give a negative result, the largest 
part of the allowed mass range of the 3z/-scheme with strong degeneracy will be excluded. 
Still a small interval (m v ~ 0.1 — 0.3 eV) will be uncovered. This will have important 
implications for the theory of the neutrino masses. 

In fig. H| we show also the upper bound on mp from data on the large scale structure of 



the Universe obtained in [[41] for values of total matter contribution to the energy density of 
the Universe, Q m = 0.4 and the reduced Hubble constant, h = 0.8. 

Without /3-decay measurements the absolute value of the neutrino mass can be deter- 
mined and the scheme can be identified provided that all the following conditions are satisfied: 

- no effect of sterile neutrinos is observed, 

- the SMA solution is established as the solution of the solar neutrino problem, 

- the neutrinoless double beta decay gives a positive result with m ee close to the present 
upper bound. 

In this case m v = m ee . However, not too much room is left for such a possibility keeping 
in mind that recent solar neutrino data disfavor the SMA solution. For all large mixing 
solutions of the solar neutrino problem, m ee gives only the lower bound on the absolute scale 
of masses. 

Let us stress that, even in the case of the SMA solution one should make the assumption 
that there are no additional contributions to the 2/30^-decay apart from the exchange of 
light Majorana neutrinos. In this scheme there are no "test equalities", that is, the rela- 
tions between m ee and the oscillation parameters, Am 2 , 9, which could allow to check this 
assumption independently. Furthermore, the determination of m ee and therefore m v will 
be restricted by uncertainties in the nuclear matrix elements. Thus, the study of the beta 
spectrum is the only way to measure the absolute scale of neutrino mass without ambiguity. 

Clearly, if the LSND result is confirmed the scheme will be excluded. 
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4 4-v schemes: Bugey, CHOOZ and LSND bounds 



Four-neutrino schemes, which explain the LSND result in terms of oscillations, have two sets 
of mass eigenstates separated by Am 2 LSND (see fig. §). Hereafter, we call them the light set 
of states and the heavy set of states. Let us consider the heavy set. The masses of states in 



this set are equal or larger than yj Am LSND . The mass differences are equal or smaller than 

Am °- or/and A < 



2 

LSND 



Both splits are much smaller than the energy resolution AE as well as masses themselves. 
So, the states in the heavy set are quasi-degenerate and their effect on the beta spectrum 
can be characterized by and p h e given in eqs. ( |I~4"D and (|20|). 



In the (2 + 2) schemes both the heavy and light sets contain two states, whereas in the 
(3 + 1) scheme one set contains 3 states, while the other set consists of only one state (see 
fig. 2). 

The v e oscillation disappearance experiments, Bugey and CHOOZ [23], impose a 
direct and very strong bound on p^, and therefore on the expected effects in /3-decay in 
all 4i/-schemes. Since Bugey and CHOOZ experiment do not resolve small mass squared 
differences, Am 2 atm and AttIq, their results can be described by 2z/-oscillations with a unique 
mass squared difference Am 2 ~ Am LSND and the effective mixing parameter 

sin 2 29 eff = 4]T \U ei \ 2 (l - J2 \Uei\ 2 ), 



where the sum runs over the heavy (or light) set. Using the definition of p e in eq. ([14]) we 
can rewrite the mixing parameter as 

sin 2 2tf e// =4p^(l-p^). (28) 

Thus, the negative results of the oscillation searches in Bugey and CHOOZ experiments give 
immediate bound on p h e as a function of Am LSND (see figs. | - )• 

For the range of masses relevant for LSND experiment (Am 2 > 0.2 eV 2 ) two possibilities 
follow from ( p8|) and the Bugey or CHOOZ bounds: 

1). Small p h e : 

p h e < 0.027. (29) 

That is, the admixture of v e in the heavy set is very small and the electron flavor is distributed 
mainly in the light set. This corresponds to the schemes with normal mass hierarchy (or to 
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normal ordering of states in the non-hierarchical schemes). Let us recall that in the schemes 
with normal hierarchy (order of states) the light set contains the pair of states which are 
separated by Am . This pair is responsible for the conversion of the solar neutrinos, and for 
brevity, we will call it "solar pair". According to (p9|), in this class of schemes one expects 
small kink at E e ~ E — nip, where vrtp > ^ /\m 2 LSND . Here, the inequality corresponds 
to the non-hierarchical case (see sect. 6). Also in the case of non-hierarchical scheme one 
predicts an observable shift of the end point associated to the masses of the light set. 
2). Large p h e : 

1 - p h e < 0.027. (30) 

The admixture of v e in the heavy set is close to one: the electron flavor is mainly distributed 
in the heavy set. This corresponds to the schemes with inverted mass hierarchy or inverted 
ordering of states. The effect in the beta spectrum consists of a "large" kink with size close 
to 1 at E e ~ E — mp. Above E e ~ E — mp the spectrum has a tail related to emission of 
the neutrinos from the light set. The rate of events in the tail is suppressed by factor given 
in eq. (|30|). If the tail is unobservable, the whole effect will look like the effect of the electron 
neutrino with a unique definite mass mp. (And it is similar to the effect in the scheme with 
three degenerate neutrinos.) 

We will further discuss these possibilities in the sections 5-7. 
Let us consider implications of the LSND result itself for the /3-decay searches. Apart from 
providing the mass scale in the range of sensitivity of future /3-decay experiment, it imposes 
an important bound on the relevant mixing parameters. 

In the 4i/-schemes under consideration the oscillations in LSND experiment are reduced 
to two neutrino oscillations with Am 2 ~ /\m 2 LSND and the effective mixing parameter 

sin 2 29 LSND = 4| £ U^l 2 = 4| £ U^U* ej \\ (31) 
j&h jel 

where summations run over the states of the heavy set in the first equality and of the light 
set in the second equality. Using Schwartz inequality we get 

sin 2 2^ D < 4(£ |t/ e ,| 2 )(£ |£/ w f) = 4p e V;, (32) 

or equivalently, 

sin 2 26 LSND < 4QT |t/«| 2 )(£ \U,j\ 2 ) = 4(1 - rf)(l - p% (33) 
iei jel 

where 

P^El^f (34) 
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is the coupling of the heavy set with the muon neutrino. The upper bound on p h ( or upper 
bound on 1 — p^, depending on the scheme) follows from CDHS experiment at high Am 2 
28fl , and from the atmospheric neutrino studies at low Am 2 [p9 . From eq. (|32l) we get a 



lower bound on p h e : 

h sin 2 2d LSND . > 

Pe> , (35) 

where both sin 2 29lsnd and are functions of Am 2 . Implications of this bound for the (3- 
decay measurements strongly depend on specific scheme, and we will discuss them in sections 
5 and 6. 

The character of the distortion of the spectrum and the sizes of the effects depend on 

1) the structure of the scheme: hierarchical or non- hierarchical; 

2) the type of hierarchy (ordering of levels): normal or inverted; 

3) the number of states in the heavy and light sets. 
In what follows we consider possible 4z/-schemes in order. 



5 Four-neutrino schemes with mass hierarchy 

In the hierarchical schemes, the masses of the states from the light set are much smaller 
than the sensitivity limit 0.3 eV (|T0|) . In the (2 + 2) schemes they are restricted by the 
atmospheric (m l < \J Am 2 tm ) or solar (m l < J Am@) neutrino mass scales. Therefore the 
observable distortion of the spectrum is only due to effect of the heavy set with the effective 
mass 

mp « ^Am 2 LSND . (36) 

As we have mentioned in the previous section, the character of distortion of the beta 
decay spectrum depends, first of all, on the type of hierarchy. 



5.1 Schemes with normal mass hierarchy 

In the schemes with the normal hierarchy (both in the (2 + 2) and (3 + 1) cases) the 
electron flavor is in the light set and p h e is strongly restricted by the Bugey result (see eq. 
(p9|) and figs. [| The beta decay spectrum has only a "small" kink with p h e < 0.027 at 
E = E - ^Am 2 LSND . 
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Additional restrictions on the /3-decay parameters may appear depending on whether the 
scheme is of the (2 + 2) or (3 + 1) type. 



1). The (3 + 1) scheme. In this scheme one gets a substantial lower bound on p h e from 
the LSND result (see eqs. (0), (|35|)). Indeed, in this scheme p h is restricted from above by 
the CDHS result . Inserting the bound on p 1 ^ into eq. (|35| ) we get the lower bound on p\ 



which is close to the upper Bugey bound or, for certain ranges of J Am 2 LSNCn even above it. 
So that only certain ranges of mp are allowed (see fig. |||). This is a manifestation of the fact 
that in the (3 + 1) scheme an explanation of the LSND result requires the v e admixture in 



the isolated state to be at the level of the upper Bugey bound ||30f| . 

Let us consider implications of the 2/50^-decay search. The contribution to m ee from the 
fourth (isolated) state dominates. It can be estimated as: 



W = jAm 2 LSND \U e4 \ 2 ~ (0.005 - 0.05) eV. (37) 



v — " 

In the hierarchical case with m 2 = J Am@, the contributions from other mass eigenstates 
can be estimated as mf) = \J ' Am 2 atm \U " e3 | 2 < 3.5- 10~ 3 eV and mfj « J 1 Am 2 & sin 2 6 Q < 7-10~ 3 
eV. Hence 

m ee « mf} = mfip h e , (38) 

Or 771/3 = "lee/Pe- 

A version of the scheme is possible in which the mass hierarchy in the light set is inverted, 



so that the states which contain the electron flavor have masses m 2 — m 3 ~ yAm 2 tm and 
mi <§C ni2- In this case we have 



m A \U ei \ 2 e lS + v/Am 2 tm (cos 2 6 Q e tQ + sin 2 9 Q ) 



(39) 



and the contribution from the light set can be comparable to the contribution from the 4th 
state. The corresponding lines are shown in fig. || According to the figure one expects 
m ee to be substantially below the present bound: We find m ee ~ 0.005 eV, ~ 0.015 eV, 
~ 0.015 -T- 0.03 eV and ~ 0.06 eV for the allowed "islands" of nip and p h e (from smallest to 
largest mp). Clearly, the observation of m ee near its present experimental bound will exclude 
the scheme. 

For the SMA solution of the solar neutrino problem, eq. (pJO) reduces to 



In principle, using this equation one can determine the relative phase S. 
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For large angle solutions of the solar neutrino problem still significant contributions can 
come from the solar pair of states and therefore two different phases (5 and a) are involved 
in the determination of m ee (see fl39D ). 



2). The (2 + 2) scheme. In the (2 + 2) scheme with normal mass hierarchy the mass 
difference of the heavy set is given by Am 2 tm , and is distributed, mostly, in the heavy set. 
According to the CHDS bound, p 1 ^ should be close to 1 in this scheme. Hence, the LSND 
bound is less restrictive: 

p h e >^sm 2 29 LSND , (40) 

and Pg can be as small as (2 — 3) • 10~ 4 (see fig. |]). The effect of the kink is unobservable 
for such a small p h e . 

Let us consider bounds from the 2/30z/-decay searches. The Majorana mass of the electron 
neutrino can be written in the following form 



( f % - I 7, ) + sin 2 (1 - \U e3 \ 2 - \U eA \ 2 ) 

mp 



(41) 



Neglecting the contribution of light neutrinos (J Am Q sin 2 9 & < 7 • 10 3 eV), we obtain from 



(fLL|) the following bound on the effective /3-decay mass: 



T. < me < (42) 

If the neutrinoless double beta decay is discovered, this inequality will put a strong lower 
bound on m^: 

m/3 > 25m ee , (43) 

where we have used the bound on p h e from the Bugey experiment. For m ee = 0.1 eV we get 
mp = 2.5 eV which is in the upper allowed region of the LSND experiment. For this reason, 
one cannot expect m ee > 0.1 eV in this scheme. 



5.2 Schemes with inverted mass hierarchy 

In the schemes with inverted hierarchy, the electron flavor is distributed mainly in the heavy 
set (see fig. §). Therefore the (3- decay spectrum should have a "large" kink with parameters 
p h e ~ 1 and mil ~ J Am 2 LSND . The light set leads to appearance of the tail at E e > E — m h n 
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with the suppressed rate determined by (1 — p£) < 0.027. To detect the signal from the tail 
one may search for the integral effect above Eq — 7714. 

Note that for the LSND region with Am 2 LSND = (6 - 8) eV 2 allowed at 99 % C.L. one 
gets m/3 > J Am 2 LSND ~ 2.5 eV which is already excluded by the Mainz result at 95 % C.L. 
(see fig. [5]). Hence, for the schemes with inverted hierarchy the preferable range of mass 
would be 

0.4 eV < nip < 1.75 eV. (44) 

1) . The (3 + 1) scheme. The upper Bugey bound and the lower LSND bound on 1 — p h e 
are shown in fig. |5|. They are similar to the bounds on p h e in schemes with normal mass 
hierarchy. 

Let us consider the implications of the 2/?0z/-decay searches. The contribution of the 
light states to m ee is negligible. As a consequence, the bounds imposed by the 2/30i/-decay 
searches are similar to the bounds in the three-neutrino scheme, and the only difference is 
that in the latter scheme, mp is a free parameter, whereas in the (3+1) scheme it equals to 
^Am 2 LSND . 

For the SMA solution of the solar neutrino problem we have m ee m J Am 2 LSND > 0.4 eV 
which is already larger than the present upper bound. So, keeping in mind the uncertainties 
of the nuclear matrix element, we can say that this possibility is disfavored. 

In fig. |5|, we show the upper bounds on mp {mp < m™ ax / 00326*0) which correspond to 
I J7"e3 1 =0 and m™ ax = 0.34 eV. The bounds are shown for different values of sin 2 2^ from 
the large mixing solution regions. Note that for the LMA solution the 99 % upper bound, 
sin 2 2# ~0.95, gives mp < 1.5 eV, which is already below the present kinematic bound. The 
best fit region sin 2 2^ Q ~ (0.6 — 0.8) leads to mp ~ (0.55 — 0.75) eV well in the range of 
the KATRIN sensitivity. For the LOW solution the maximal mixing is possible for which 
008 26*0 = 0, and the whole region of mp up to the present kinematic bound @ is accepted. 

2) . The (2 + 2) scheme. Here the effects are similar to those in the (3 + 1) scheme with 
two differences: 

i) . As it is shown in fig. || the LSND result is less restrictive and the rate in the tail can 
be substantially lower. 

ii) . The bound on mp from the neutrinoless double beta decay coincides with the bound 
in the (3 + 1) scheme at \U e ^\ = 0: 

^ee / . _ \ 

m ee <mp < — . (45) 

cos 20 Q 



17 



6 Non-hierarchical schemes 



Let us consider schemes in which the masses of states from the light set are also in the 
range of sensitivity of KATRIN experiment: mi > 0.3 eV. Clearly, these states are quasi- 
degenerate, and their effect on the /3-decay spectrum can be characterized by the effective 
mass, Trig, and the total coupling with the electron neutrino, p e . From the unitarity condition 
we have p l e + p\ = 1. 

In the non-hierarchical schemes, the effect of neutrino mass on the /3-decay spectrum 
consists of the kink and the shift of the end point. The kink is at 

E kink = Eq _ m h^ ( 46 ) 

where 

mj = y/Aml SND + (m l p¥, (47) 

and its size is determined by p h e . In contrast to the hierarchical case, the mass of the heavy 
set is not fixed by \j 'Amj jSND . The tail above the kink, E e > E* mk , is described by p\ and 
the end point is shifted to 

E? = Eo- m l g . (48) 

Note that if m l g ~ 0.5 eV, from ([IT]) we obtain the effective mass of the heavy set 
mg ~ 2.5, 1.5 and 0.8 eV for Am 2 LSND = 6, 1.5 and 0.4 eV 2 correspondingly. For large 
values of /S.m 2 LSND both the structures (the kink and the bending of the spectrum at the 
end point) can in principle be separately detected by KATRIN. For /\m 2 LSND at the lower 
allowed end the difference of the two effective masses becomes small: rrig — m l g < 0.3 eV, 
and these two structures may not be resolved. 

With increase of m l g the scheme transforms into the scheme with four degenerate neutri- 
nos. Already at m l g = 1 eV, we get the mass of the heavy set m h g = 2.7, 1.5, and 1.2 eV, 
for /S.m 2 LSND = 6, 1.5, and 0.4 eV 2 . For the smallest mass, m^=1.2 eV, the difference of 
masses, m!g — m l B < 0.2 eV, is too small to be resolved. With the increase of m l g the two 
structures in the spectrum, the kink and the bending, merge. 

The parameters of the kink and of the tail depend on specific properties of the scheme. 
We will call it the scheme with normal ordering of states when the electron neutrino is 
distributed mainly in the light set. And we will refer to the opposite situation, when v e is 
in the heavy set, as to the scheme with inverted order of states. 

In the non-hierarchical schemes the mass spectrum is shifted to larger values of mass. 
The oscillation pattern, however, is not changed and the oscillation bounds are the same as 
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in the hierarchical schemes described above (see fig. |3], |5|). However, the implications of the 
2/50z/-decay searches are changed. 

6.1 Schemes with normal order of states 

The electron neutrino is mainly in the light set. The beta decay spectrum should have a 
small kink at Eq — with the size restricted by Bugey experiment: p h e < 0.027, and a 
strong bending of spectrum at E — m^. 

In the fig. |7| we show the bounds on the /3-decay parameters and p h e for two represen- 
tative values of the mass squared difference: Am 2 LSND = 1.75 eV 2 , which corresponds to the 
weakest bound on p h e from Bugey experiment, and Am 2 LSND = 0.227 eV 2 from the lowest (in 
Am 2 scale) LSND region. 

For Am 2 LSND = 1.75 eV 2 the allowed region of the /3-decay parameters is between the 
vertical lines at p h e = 0.013 (dashed) and p h e = 0.026 (solid) for the (3 + 1) scheme (shad- 
owed), while for the (2 + 2) scheme the valid region stays between p h e = 0.004 (dash-dotted) 
and p h e = 0.026 (solid). 

For Am 2 LSND = 0.227 eV 2 the allowed regions are substantially smaller: for the (3 + 1) 
scheme the region is between lines at p^ = 0.0095 (dashed) and 0.010 (solid). For the (3 
+ 1) scheme the region (shadowed) is restricted by lines at 0.0002 (dash-dotted) and 0.010 
(solid). All the regions are bounded from above by the Mainz result. 

The implications of the 2/30^-decay searches depend on the specific arrangements of 
levels. 

1. (3 + 1) scheme. The contribution to the effective Majorana mass, m ee , from different 
mass eigenstates can be evaluated in the following way. The solar pair of states (u e is mainly 
in these states) yields the contribution: 

™ { ee +2) » m™ n = m l p {\ - |f/ e3 | 2 - |t/ e4 r)(cos 2 6 Q + e lS sin 2 ) . (49) 

The contributions from the two other states are 

mg> = m l pU 2 e3 , mfj = mJ[/ 2 4 . (50) 

In general these contributions may have arbitrary relative phases and cancel each other in 
the sum. The contribution from the third level, which belongs to the light set, is restricted 
by the CHOOZ bound: mf) < 0.05 eV for m 1 ^ < 1 eV. In turn, the fourth contribution (from 
the isolated level) is restricted by the Bugey result: rrvg < 0.05 eV. The contribution of the 
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solar pair ( |4"9"| ) depends on the solution of the solar neutrino problem and, in most of the 
cases, dominates over other contributions. Indeed, for the SMA solution we get m™" « m 1 ^ > 
0.3 eV (if masses of the light set are in the range of sensitivity of KATRIN experiment). 
For the LMA solution the cancellation of the two terms in (0^) may occur but typically, 
m™" ~ (0.2 — l)mi > 0.1 eV is large enough to be detected in the forthcoming 2/30i/-decay 
experiments. For other solutions with large mixing angle (LOW, VO) the cancellation can 
be stronger. The solar pair contribution can be comparable with two others if the mixing is 
close to maximal: sin 2 29 & > 0.98 and the two terms in ( fi-9|) have opposite signs. 

So, in general we expect the effective mass of the Majorana neutrino to be m ee ~ 0.1 eV, 
that is, not too far from the present experimental bound. The identification of the solution 
of the solar neutrino problem can clarify a situation leading to more definite predictions. 

If the solar pair gives the dominant contribution, the measurements of m ee , m 1 ^ and 
6 Q will allow to determine the relative phase of masses in the solar pair, 5 (see eq. (PS|)). 
Otherwise, due to the presence of three different phases we cannot determine values of these 
phases, and m ee will give only a lower bound on the mass scale. 

Assuming the maximal cancellation of contributions in m ee we find from (f49|) and (|50|). 
an implicit upper bound on mjg as a function of p h e : 



m. 



p h e )\cos26 Q \ - J(m l y + Aml SND p h e - mUl - | cos2# |)|[/ e3 | 2 



(51) 



We show this bound on mp for different values of sin 26 Q in the fig. [7|. Note that for 
sin 2 ^© >0.95, the effect of non-zero \U e 3\ 2 (the last term in the left hand side of eq. (|5I])) is 
nonnegligible but it decreases with sin 2 26 l 0. In fig. [7|, the line marked by "3+1" shows the 



2/?0z/-decay bound for the (3 + 1) scheme at |?7 e3 | 2 =0.05 and \J Am 2 LSND =0.477 eV. The 
lower line in the pair marked by "2+2" is the corresponding bound for the (3 + 1) scheme 
at \U e 3\ 2 = 0. At \U e 3\ 2 = the bounds for the (2 + 2) and (3 + 1) schemes coincide (see 
below). The pair of lines marked by "2+2" illustrates dependence of bound on Am 2 LSND . 



The upper line (a) is for J Am 2 LSND = 1.32 eV while the lower one is for 0.477 eV. For 



smaller values of sin 2 2# , the lines have been calculated at |?7 e 3| 2 =0 and y Am 2 SAri) =1.32 
eV. 

2. (2 + 2) scheme. Now the heavy set contains two states, u 3 and z/ 4 , and the solar 
pair is in the light set. The contribution to the Majorana mass, m ee , from the light set is 
described by eq. (|49|). The heavy set gives 



mg +4) = < 



U e z\ 2 + e iS \U e ,\ 



< ™y e (52) 
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which is restricted by the Bugey bound: < 0.04 eV for wJl < 0.1 eV. Again, the solar 

pair gives the dominating contribution unless the solar mixing is very close to maximal and 
the phase in ( |49|) is close to n. 

Assuming the maximal cancellation of contributions in m ee we find from fl52|) and (f49[), 
an implicit upper bound on as a function of p^: 



mj,(l - p h e )\cos26 Q \ - ^(m^ + Aml SND p h e 



< m ee . (53) 



We show this bound on mp for different values of sin 2 29 in the fig. |7|. The pair of lines 



marked by "2+2", corresponds to two different values of J Am 2 LSND : 1.32 eV (upper line) 



and 0.477 eV (lower line). For other values of sin 2 2^ , J Am 2 LSND is taken to be 1.32 eV. 
Note that for the (3+1) scheme with \U e3 \ 2 = 0, the bounds from the 2/30z/-decay searches 
are the same as for the (2+2) scheme. 



6.2 Schemes with inverted order of states 

The electron neutrino is mainly distributed in the heavy set, so that one expects a large kink 
at Eq — rrip with the size close to 1 and the suppressed tail with the end point at Eq — mi. 
In the tail the rate is restricted by the Bugey bound: p l e < 0.027. As for the case of the 
spectrum in the schemes with inverted mass hierarchy, we can conclude that Am 2 LSND < 3 
eV 2 . 

In fig. |8|, we show the bounds on the relevant beta decay parameters: m^, the effective 
mass of the heavy set, and (1 — p^) which determines small admixture of the electron neutrino 
in the light set. Clearly, 

™ 1 P = ^(™p) 2 -&™LSND (54) 

and 

m h p > ^Am 2 LSND . (55) 

We show the bounds for two representative values of \J /S.m 2 LSND : 1.32 and 0.447 eV (as 
for the scheme with normal ordering). The allowed ranges for (1 — p^) determined by the 
oscillation experiments are the same as the ranges of p h e for normal ordering (fig. ||). The 
allowed regions for the (3 + 1) scheme are shadowed. The allowed values of are restricted 
by the Mainz limit from above and by the LSND result from below (see eq. (|55|)). 

Let us consider the bounds from the 2/30i/-decay searches. The effective Majorana mass 
can be immediately obtained from the results for the scheme with normal ordering by the 
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interchange: «-> m^. This leads to a stronger dominance of the solar pair contribution 
to the effective Majorana mass. 

1. (2 + 2) scheme. The contribution from the heavy set is given by 

m h ee ~ mZ n = m h p{± - \U e3 \ 2 - |*7 e4 | 2 )(cos 2 6 Q + e tS sin 2 ), (56) 
whereas the contribution from the light states can be written as 

m l ee = m l 3 \U 2 el + e lS 'U 2 e2 \ < m l g p l e , (57) 

and according to the Bugey result: p\ < 0.027. Implications of the 2/30z/ searches are similar 
to those in the hierarchical case. However now m^, and consequently m ee , can be even larger 
than JAm 2 LSND . 



From eqs. fl56] , |57 ) we find the lower bound on m ee : 

m ee > m h p{l - p l e )\ cos26'0| - m l pp l e . (58) 

Using this inequality and the upper experimental bound on m ee , we get an implicit upper 
bound on wig as a function of p h e : 



m y e \ cos2# | - ^{m h 3 y - Aml SND (l - p h e ) < = 0.34 eV. (59) 

The bounds for different values of the solar mixing parameter are shown in fig. || The 
identification of the solution of the solar neutrino problem and measurements of sin 2 29 Q as 
well as mild improvement of the bound on the Majorana mass will have strong impact on 
this scheme. For instance, as follows from the fig. || the possible bounds: sin 2 29 & < 0.9 
and m ee < 0.1 eV would exclude whole the region of parameters of the scheme down to the 
KATRIN sensitivity limit. 

2. (3 + 1) scheme. The contributions to m ee from the heavy and the light sets are equal 

to: 



m h ee = m h p 



1 - \U e3 \ 2 - \U e4 \ 2 )(cos 2 6 Q + e td sin 2 ) + U 2 3 (60) 



and 

m l ee = m l 3 (l-p h e ), (61) 

respectively. In ( |50D \U e3 \ 2 is restricted by the CHOOZ results and m l ee is restricted by the 
Bugey results: m l ee ~ 0.03 eV taking m l B < 1 eV. The contribution from the heavy set is 
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similar to the one in the scheme with three degenerate neutrinos. For the largest part of the 
allowed parameter space the contribution of the solar pair to m ee dominates. Still significant 
cancellation is not excluded which can cause the two contributions in (|60f) (from the solar 
pair and the third mass eigenstate) to be comparable. Note that when \U e 3\ 2 is smaller than 
0.015 — 0.05, m l ee can be as large as m^\U e3 \ 2 . 

For the SMA solution we have m ee > 0.6 eV for mi > 0.3 eV, so that such a possbility is 
excluded by the present bound (||). 

Assuming the maximal cancellation of contributions in m ee , we find from 
an implicit upper bound on as a function of p h e : 

m h p p h e \ cos20 | - ^/(m^) 2 - Am 2 LSND (l - p h e ) - mj(l - | cos2£ |)|£/ e3 | 2 < m ee . (62) 

We show this bound on for different values of sin 2 29 & in the fig. |8|. Note that for 
sin 2 200=0.95, the effect of non-zero \U e3 \ 2 is non-negligible but for smaller values of sin 2 29 Q 
we can neglect |t/ e 3| 2 . In the fig. ||, for sin 2 200=0.95, we have taken |[/ e 3| 2 =0.05 and for 
other values of sin 2 26*0 we have set \U e 3\ 2 equal to zero. 



6.3 4z/- schemes without LSND 

Apart from the LSND result, there is a number of other motivations to introduce new 
neutrino mass eigenstates. In particular, the sterile neutrino in the eV-range has been 
discussed in connection to the supernova nucleosynthesis (r-processes) |j31| . The mixing of 
the keV-mass sterile neutrino with the active neutrinos can provide a mechanism of the 



pulsar kicks ||32| . The keV sterile neutrinos may compose the warm dark matter of the 
Universe ||33|| . Small mixing of the sterile neutrino with the active neutrinos can induce the 
large mixing among the active neutrinos . 

Light (SU(2) x U(l)) singlet fermions ("sterile neutrinos") can originate from some new 
sectors of the theory beyond the standard model. Neutrinos, due to their neutrality are 
unique particles which can mix with these fermions. So, searches for the effects of the sterile 
neutrinos are of fundamental importance even if these fermions do not solve directly any 
known problem and thus their existence is not explicitly motivated. 

In this connection we will consider a general four neutrino scheme in which three (dom- 
inantly active) neutrinos are light and the fourth neutrino (dominantly sterile) has a mass 
in the eV - keV range. The three light neutrinos may form a hierarchical structure with the 
heaviest component being below 0.07 eV, so that their masses will not show up in the plan- 
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ning beta decay experiments. The fourth neutrino has a small mixing with active neutrinos, 
and in particular, with the electron neutrino. 

The scheme is similar to the (3 + 1) schemes with normal mass hierarchy. The difference 
is that now the mass m 4 and the mass squared differences m| — mf (i = 1,2,3) are not 
restricted by the LSND result so that m 4 can be larger or much larger than 2.5 eV. 

Let us consider the possible effect of this fourth neutrino in the beta decay. We con- 
centrate on the range of masses ~ (0.5 - 5) eV which satisfy the cosmological bounds. The 
fourth state produces the kink in the beta decay spectrum at E = E — m 4 with the size 

Pe = \U e4 \ 2 . (63) 

Let us evaluate the allowed range of p e for different values of m 4 . At m 4 > 1.5 eV the 
strongest bound follows from the CHOOZ result: p e < 0.027 (we assume that other neutrinos 
are much lighter, but generalization to the non-hierarchical case is straightforward). This 
bound does not depend on the mass for m > 1.6 eV. 

Another direct restriction comes from the 2/50i/-decay searches provided that neutrino is 
the Majorana particle: 

m 4 < (64) 

Pe 

where m™ ax is the upper bound on the effective Majorana mass of the electron neutrino. 
Taking the present bound m™ ax =0.34 eV and the maximal allowed value of p e we find that 
m 4 < 12.6 eV. Thus, one may see the kink of the 3% size (or smaller) in the energy interval 
(0 — 13) eV below the end point. The recent cosmological bound ([I]) shrinks substantially 
this interval. 

Additional restrictions on the possible effects appear if there is a substantial admixture 
of the muon flavor in the fourth state. In this case one predicts the existence of the — v e 
oscillations with the effective mixing parameter 

sin 2 2# eM = 4|?7 e4 | 2 |[/ M | 2 (65) 

and Am 2 « m\ = (1-100) eV 2 . For Am 2 > 7 eV 2 the stronger bound, sin 2 29 ei _ L < 1.3 -10~ 3 , 
is given by KARMEN experiment ||35|| . Clearly, these bounds are satisfied, if |t^4| 2 is small 
enough. However, if (t^l 2 > \U e ^\ 2 (which might be rather natural assumption) the bound 
from KARMEN experiment becomes important. Taking l^^l 2 = |t/ e4 | 2 we get from (|65| ) 

p e < ^y/sin 2 20 e „, (66) 

and for m 2 > 7 eV 2 , it follows from (|66|) that p e < 1.7 • 10 -2 which is stronger than the 
Bugey bound. For heavy neutrinos (in keV range) the neutrinoless double beta decay gives 
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very strong bound on the size of the kink: p e < 3 • 10 _4 (lkeV/ni4) which will be very difficult 
to observe. This bound does not exist if the keV neutrino is the Dirac (or pseudo-Dirac) 
particle. 

7 Beta decay measurements and the neutrino mass 
spectrum 

In this section we consider the astrophysical and cosmological bounds on the neutrino 
mass. We will discuss possible future developments in the field. 

7.1 Beta decay measurements and supernova neutrinos 

Studies of the supernova neutrinos open unique possibility to test the schemes of neutrino 



mass and mixing ||45|| . Therefore they may have an important impact on predictions for future 
beta decay measurements. 

Considering the level-crossing patterns j|5] and the adiabaticity conditions in various 
resonances it is easy to show that in all the (2 + 2) schemes with inverted mass hierarchy 
(or ordering of the states) the originally produced i/ e -flux is almost completely converted to 
some combination of v T and ^-fluxes at high densities in the resonance associated to 
Am 2 LSND . The mixing parameter in this resonance, given by sin 2 2Q LS ndi is large enough to 
garantee the adiabaticity of the conversion. In this case the v e survival probability equals to 
P m sin 

Olsnd < 10~ 2 . At the same time, the z/ e -flux observed at the Earth appears as a 
result of conversion 

i/p, v T -> v e (67) 

at high densities inside the star. Therefore the spectrum at the Earth will practically coincide 
with the hard original spectrum of and v T : 

F- e {E) « (68) 

and moreover, this result does not depend on the solution of the solar neutrino problem. Such 
a hard spectrum of v e is strongly disfavored by the SN1987A data pEEfl . Future detections of 
the Galactic supernovae can exclude the conversion (|67|) , and consequently the schemes will 
be excluded, completely. 

In the (3+1) scheme with inverted mass hierarchy (ordering) the result of conversion 
depends on the solution of the solar neutrino problem [001. As in the (2+2) scheme, the 
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original z/ e -flux is converted to a combination of v T and £/ s -fluxes. For the SMA solution no 
opposite conversion (that is, v T and to u e ) occurs. Therefore, in this scheme the tran- 
sitions lead to practically complete disappearance of the z/ e -flux. The suppression factor, 
sin 2 Olsnd < 10~ 2 , can not be compensated by the allowed increase of the original flux. The 
disappearance of u e contradicts the data from SN1987A, so that the scheme is excluded. No- 
tice that this scheme is also practically excluded by the present bound from the neutrinoless 
double beta decay. 

In the case of the LMA solution some part of the original v^— and v T — fluxes will be 
transformed to the v e — flux at low densities, so that at the surface of the Earth one expects: 

F S (E) ps sin 2 8 Q Fp(E). (69) 

Thus, u e will have hard spectrum suppressed by factor 1/3 — 1/2. This is again disfavored 
by the SN1987A data. 

Notice that in all these schemes the electron neutrino flux is not elliminated from the 



region proposed for the r-processes |31| . So that the mechanism of production of the heavy 
elements will not work. 

In contrast, the non-hierarchical 4z/-schemes with normal ordering are well consistent 
with the SN1987A data and they predict an observable effect in the /3-decay spectrum, as 
was discussed in sect. 6. 

7.2 Beta decay and forthcoming experiments 

The results of the forthcoming oscillation as well as non-oscillation experiments can substan- 
tially influence the both predictions of the effects of neutrino mass and mixing in the beta 
decay spectrum and the significance of future beta decay measurements. In particular, 

(i) the identification of the solution of the ^ -problem and measurements of relevant 
oscillation parameters, 

(ii) the MiniBooNE result, 

(iii) further searches for the neutrinoless double beta decay 

will have crucial impact. Also further improvements of the bound on |i7 e 3| will be important. 
Cosmology can give a hint for the absolute scale of neutrino mass. 

Let us analyze consequences of possible results from these experiments. 

1). The solution of the solar neutrino problem can be identified in the forthcoming 



experiments: SNO |§, g7]], KamLAND |8[, BOREXINO ||. The identification will not 



influence the predictions for the /3-decay immediately. Indeed, from these experiments we 
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will get specific values (ranges) of A?tIq and mixing angle 9 & , i.e., the distribution of the 
electron flavor in the solar pair of states will be determined. But, the effects in the /3-decay 
are not sensitive to a particular value of Amg, since for all possible solutions Am|, can 
not be resolved in /3-decay searches. Also the effects are not sensitive to the distribution 
of the electron flavor since they are determined by the sum over states in the solar pair: 
\U e i\ 2 + \U e 2\ 2 ~ 1- However, the identification of the solution of the solar neutrino problem 
will influence substantially the bounds on the /3-decay effects from the 2/30^-decay searches. 
A number of schemes discussed here will be excluded and for other schemes the possible 
effects in the /3-decay spectrum will be strongly restricted. 

The key issues are whether the correct solution of the solar neutrino problem is the small 
mixing solution or the large mixing solution, and if it is the large mixing how large the 
deviation from maximal mixing is. 

Suppose that the SMA solution will be identified, then the following information can 
be obtained from the 2/30i/-decay searches in the assumption that the Majorana neutrino 
exchange is the only mechanism of the decay: 

• For the 3^-scheme this will imply that mp < 0.34 eV and further moderate improve- 
ment of the 2/30^-decay bound will exclude the scheme. 

• According to present data: J Amj jSND > 0.38 eV (at 99 % C.L.) |Hfl . Therefore in 
the schemes with inverted mass hierarchy or inverted order of states we get m ee > 
J Am| 5JV£) > 0.38 eV. On the other hand, the bound from the 2/30z/-decay is m ee < 
0.34 eV (at 90 % C.L.). Therefore these schemes are excluded at stronger than 90 
% C.L. One should, however, keep in mind the uncertainties of the nuclear matrix 
elements. Future double beta decay measurements will be able to confirm and improve 
the bound. Similar conclusion can be made for non-hierarchical schemes with normal 
ordering of states and mi > 0.34 eV. 

Thus, the schemes which will survive after the identification of the SMA solution, are 
the schemes with normal hierarchy as well as the non-hierarchical schemes with normal 
order of states and mi < 0.34 eV. 

• If m ee turns out to be close to the present bound (e.g., ~ 0.2 eV) and the LSND result is 
confirmed, the only possibility will be the non-hierarchical scheme with normal ordering 
of states and mi = m ee . 

Suppose now that one of the large mixing solutions of the solar neutrino problem will 
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be identified. In this case a possible cancellation between various contributions to m ee will 
relax the bounds from the 2/?0^-decay. 

• The important point is that the present upper bound on m ee is already smaller than 
^Am 2 LSND : 

™? e ax < ^mlsND 

(although, one should keep in mind the uncertainties of the nuclear matrix element). 
At the same time, > J Am 2 LSND . Consequently, the schemes with inverted mass 
hierarchy or inverted ordering of states require cancellation between the contributions 
from the solar pair states, and therefore, there should be CP-violating phase difference 
in the corresponding mass eigenvalues. The scheme with the SMA solution of the solar 
neutrino problem is practically excluded. 

The upper bound on sin 2 29 & (lower bound on the deviation from maximal mixing) 
restricts a possible cancellation of contributions to m ee . This, in turn, gives an upper 
bound on the mass m^. Therefore, further improvements of the bounds on m ee and 
sin 2 29q can strongly restrict the parameter space of the schemes and even exclude 
them. 

Similar conclusions hold for the non-hierarchical schemes with m 1 ^ > m™ x ~ 0.34 eV. 

• As we discussed in section 5, the Majorana mass m ee gives the lower bound on the 
mass of the solar pair. So, if m ee turns out to be close to the present bound (f|), we 
can exclude the schemes with normal hierarchy (in which m ee < m l < 0.07 eV). 

• If m ee is much smaller than the present bound (e.g., ~ 0.01 eV) the scheme should 
have the normal mass hierarchy or should lead to a strong cancellation of contributions 
to the 2/90^-decay. 

2). MiniBooNE experiment will give strong discrimination among the possibilities. 

If MiniBooNE does not confirm the LSND result, a large class of schemes discussed 
here will be excluded. There will be no strong motivation to consider the Av schemes (see, 
however, sec. 6.3). 

Also further improvements of bounds on the involvement of a sterile neutrino in the 
solar and in the atmospheric neutrino conversions may give independent confirmation of the 
3z/-schemes. 
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We will be left with the three neutrino scheme with strong degeneracy or with schemes 
having more than three mass eigenstates without observable signal in MiniBooNE (see sect. 



So, in this case the searches of the neutrino mass in the sub-eV range will basically test 
the Zv scheme with strong degeneracy. The observation of the shift of end point to E — mp 
in KATRIN experiment will be the proof of the scheme with mi ~ ni2 — ~ mp. 

Further insight can be obtained confronting the results of the /?-decay and the 2/50i/ decay 
measurements, as we have discussed in sect. 3. 

If MiniBooNE confirms the LSND result we will be forced to consider the 4z/-schemes. 
In this case we get for the absolute value of the mass of heavy set: 



where the equality corresponds to the scheme with the mass hierarchy {wJl ^> ml). 

MiniBooNE experiment will not only check the LSND result but also further restrict the 
oscillation parameters. Moreover, it may allow to disentangle the (3+1) and (2+2) schemes. 
Further searches of the sterile neutrinos in the solar and atmospheric neutrino fluxes should 
discriminate the (2+2) and (3+1) schemes ||30|| . 

3). Let us consider implications of future cosmological measurements. The present and 
future cosmological bounds on the mass scale of neutrinos are summarized in the table [1| 
The bound on m u taken from |4l| corresponds to the energy density of matter fl m = 0.4 and 
the reduced Hubble constant h = 0.8. 

Note that by chance the best cosmological bound (]l|) coincides numerically with the best 
laboratory limit (|8J). However, in contrast to the latter, the cosmological bound is valid for 
any flavor including the sterile neutrino, provided that this neutrino had been equilibrated 
in the Early Universe. 

The bound (jl]) (obtained for one neutrino in the eV range) applies immediately to the 
(3 + 1) scheme with normal hierarchy. In this scheme the heaviest (isolated) state can 
produce only a small kink in the /3-decay spectrum which will be difficult to detect. So the 
cosmological bound being confronted with the value of J /S.m 2 LSND will play important role 
in checking the scheme. The bound is even stronger for the non-hierarchical (3 + 1) schemes. 

The bound is also immediately applied to the schemes with additional heavy neutrino 
(sect. 6.3). Even for very small admixture of the electron neutrino in this state (|t/ e 4| 2 <C 
0.05) this, predominantly sterile, neutrino will have equilibrium concentration in the Uni- 
verse. 



6). 




(70) 
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Table 1: Cosmological bounds on the neutrino mass scale. All the four- neutrino schemes are 
considered to be hierarchical. 



Mass scheme 


present m v eV, ref. 


future m u eV, ref. 


3 v or (3+1) 
inverted 


1.8 f 
2.5 


11 

14 


1 


0.41 


12 




(2+2) 
inverted or normal 


3.0 
3.8 


13 
14 


1 
1 


0.57 






(3+1) normal 


5.5 [ 

7.6 [ 
2.2 


13 
14 

1 


J 

1 


0.99 







The bound on mass of two or three heavy degenerate neutrinos will be stronger than 
(|I|). However, the decrease of the limit is weaker than just n~ l , where n is the number of 
the degenerate neutrinos. We estimate the bounds for two and three neutrinos performing 
rescaling of the bound (|1]) according to results for 1, 2, and 3 neutrinos in |[42|| . Thus, for 
two degenerate neutrinos we get from [43j and | m < 1.3 eV. In the (2 + 2) schemes with 
normal mass hierarchy or with normal ordering this limit excludes the upper "island" of 
parameters allowed by LSND (see fig. f|). In the case of inverted hierarchy it confirms the 
Mainz result. 

For the three degenerate neutrinos we find m < 0.9 eV. This bound excludes significant 
parts of the otherwise allowed regions of the 3i/-scheme and the (3 + 1) schemes with inverted 
hierarchy or order of the states. 

Forthcoming cosmological data can further substantially improve the bounds. In the last 
column of the table 1 we show the bounds which can be obtained using data from the Sloan 
Digital Sky Survey (SDSS) @ for Vt m h 2 < 0.17. 



7.3 Beta decay measurements and the neutrino mass spectrum 

The identification power of the /3-decay studies depends on whether future measurements 
will be able to observe the small kink and the tail after large kink or not. 

If the /3-decay measurements are not able to identify the "small" kinks and the suppressed 
tail, the only expected distortion effect is a shift of the end point and the corresponding 
bending of the spectrum. 

Suppose that the LSND result will not be confirmed and the effects of sterile neutrinos 
will not be not found neither in the solar, nor in the atmospheric neutrino fluxes. This will 
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be the evidence for the three neutrino scenario. (The existence of additional sterile neutrino 
which mixes weakly with the block of active neutrinos does not change the conclusion.) As 
we have discussed, in this case the /3-decay parameter, mp, will determine the scale of three 
degenerate neutrinos. 

If the LSND result is confirmed and Am 2 LSND is measured, the important conclusions 



will be drawn from the comparison of the values of mp and y Am 2 LSND . 

Let us remind that experiments which are insensitive to the small kinks and tails will find 



mp = for the schemes with normal mass hierarchy, and mp = J Am 2 LSND for the schemes 
with inverted mass hierarchy. Any type of relation is possible for the non-hierarchical scheme 



with normal ordering of states: (i) mp < J Am 2 LSND , (ii) mp > J Am 2 LSND in the case of 



relatively small Am 2 LSND or (iii) mp « J Am 2 LSND (which looks as an accidental coincidence 
since this will imply the equality m 1 = J{m h ) 2 — [m 1 ) 2 ). For the non-hierarchical schemes 



with inverted order one has mp > \J Am 2 LSND . 
Therefore: 

• A negative result of the neutrino mass measurements will be the evidence of the scheme 
with normal hierarchy or the non-hierarchical scheme with normal order and m l below 



the sensitivity limit: m < 0.3 eV (see eq. (10)) 



If it is established that mp < ^ Am 2 LSND , the non-hierarchical scheme with normal 
order of states should be selected. Moreover, the mass scales will be completely deter- 
mined: m l = mp, and m h = Jm 2 + Am 2 LSND . 



The inequality mp > Am LSND will testify for the non-hierarchical scheme with 
inverted order of states. In this case: m h = mp and m l = \Jm 2 p — Am 2 LSND . The 
inequality can correspond also to the non-hierarchical scheme with normal order, so 



that mp = m l and m h = J ml + Am 2 LSND 



If mp coincides (within the error bars) with Am 2 LSND , one of the following possibilities 
will be realized : 



a) the schemes with inverted mass hierarchy and mp ~ m^ ~ y Am 2 LSND , 

b) the non-hierarchical scheme with inverted order of levels and relatively small mass 
of the light set, m l , so that m h is only slightly (within the error bars) larger than 



/ ^ m LSND- 



c) non-hierarchical schemes with normal order of levels and mp —m. In this case the 



equality m l ~ J Am 2 LSND is accidental. 
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Recently, possible implications of results from LSND and KATRIN experiments have been 
discussed also in fj?fl . 



If future /3-decay experiments detect the "small" kinks and the suppressed tail, we will be 
able to measure p h e and p l e and unambiguously discriminate the hierarchical schemes from the 
non-hierarchical ones and the inverted scenarios from the normal ones even without using 
the LSND result. This will be independent test of the 4z/-scheme. Indeed: 

• The spectrum with the small kink at nip and the tail without bending, will be the 
evidence of the scheme with normal mass hierarchy. 

• The spectrum with the large kink at mp and the tail without bending will testify for 
the schemes with inverted mass hierarchy. 

• The spectrum with the small kink and the bending at higher energies will correspond 
to the non-hierarchical scheme with normal order of the states. 

• The spectrum with the large kink and the suppressed tail with shifted end point (the 
latter will probably be impossible to establish) will indicate the non-hierarchical scheme 
with inverted order of states. 

Notice that all these results are the same for the (3+1) and (2+2) schemes. The study of 
the /3-decay spectrum cannot distinguish these schemes. The (3+1) and (2+2) schemes can 
be discriminated by studies of effects of the sterile neutrinos in the solar and atmospheric 
neutrino oscillations, by MiniBooNE experiment, by searches for the oscillations in the v e 
and disappearance channels, etc. 

The detection of the small kinks will also open a possibility to search for the mixing of 
the sterile neutrinos which are not associated with the LSND result. 



7.4 Measuring the absolute mass scale 

Without direct kinematic measurements, the absolute scale of neutrino mass can be es- 
tablished only in certain exceptional cases. This will be possible if solar neutrino data are 
explained by the SMA solution and the 2/30z/-decay is observed. In this case m ee will give the 
mass of the solar pair: m sun = m ee . Then using the oscillation results one can reconstruct 
the whole spectrum. 

In the 3z/-scheme m ee will immediately determine the mass of all the three quasi-degenerate 
neutrinos. 
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In 4z/-schemes the mass reconstruction will depend on the type of the scheme: 

In the scheme with normal mass hierarchy m ee is expected to be small: m ee < \J Am^ m < 

0. 07 eV. So, if a small m ee is detected or a strong bound on m ee is obtained, the mass of 
the heavy set will be given by the LSND result: m h yj Am 2 LSND . Establishing the scale 
of light masses still will be rather problematic. Moreover, if the 2/30z/-decay searches give 
only the upper bound on m ee , we should assume that neutrinos are the Majorana particles 
to make conclusion on the mass. 

In the non-hierarchical scheme with normal ordering of states, m ee can be as large as the 
present upper bound. Thus, m ee will determine the mass of the light set: m l = m ee , and for 
the heavy set we find m h = \jm\ e + Am 2 LSND . 

The schemes with inverted mass hierarchy or inverted ordering are almost excluded by 
the fact that already present data indicate inequality (m™ ax ) 2 < Am 2 LSND . 

For the large mixing solutions of the solar neutrino problem (LMA, LOW, VAC) the 
absolute mass scale can not be restored from m ee due to possible cancellation which depends 
on unknown CP-violating phase. Inversely, the data from the /5-decay measurements can be 
used to determine this phase. If m ee is measured, we will be able to put both the lower and 
the upper bounds on the absolute scale of masses: mp < m ee / cos2# and mp > m ee . 

Even in those special cases where the determination of the absolute scale is possible there 
are two problems: 

• Uncertainties of the nuclear matrix elements will lead to significant uncertainty in the 
determination of the absolute scale. 

• We should assume that the exchange of the light Majorana neutrinos is the only mech- 
anism of the 2/30z/-decay. 

In view of this, and keeping also in mind that the SMA solution is disfavored by the 
latest solar neutrino data, we can conclude that developments of the direct methods of 
determination of the neutrino mass (and KATRIN may be only the first step) is unavoidable 
if we want to reconstruct neutrino mass spectrum completely. 

8 Conclusions 

1. We have studied effects of the neutrino mass and mixing on the /3-decay spectrum in 
three neutrino schemes which explain the solar and atmospheric neutrino data as well as in 
all possible 4z/-schemes which explain also the LSND result. 
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We find that schemes which can produce an observable effect in the planned sub-eV 
measurements should contain the sets (one or more) of quasi-degenerate states. The only 
exception is the (3 + 1) scheme with normal mass hierarchy. However it leads to a small 
kink which will be difficult to observe. We show that the effects in the /3-decay spectrum 
are described by the effective masses of the quasi-degenerate sets, and their coupling 

with the electron neutrino, 

2. At present, a rather wide class of realistic schemes exist which can lead to an observable 
effect in the sub-eV studies of the /5-decay spectrum. We show however that future oscillation 
experiments and 2/90z/-decay searches can significantly restrict these possibilities. 

3. The three neutrino schemes which explain the solar and atmospheric neutrino anoma- 
lies in terms of neutrino oscillations can lead to an observable effect in future /3-decay mea- 
surements only if all three mass eigenstates are quasi-degenerate. The /5-decay measurements 
give unique possibility to identify these schemes. Even if SMA solution is established and 
the neutrinoless double beta decay is discovered, so that m ee sets the scale of the neutrino 
masses, the question will rise whether the neutrino exchange is the only mechanism of the 
neutrinoless double beta decay. Only comparison of the m ee with results of the beta decay 
measurements will give the answer. In the case of large mixing solutions of the solar neutrino 
problem, simultaneous measurements of the m ee and mp as well as the solar mixing angle 
open the possibility to determine the relative CP-violating phase of the mass eigenvalues. 

4. In the four neutrino schemes which explain also the LSND result, observable effects 
in the /5-decay are strongly restricted by the Bugey and CHOOZ bounds. Four types of 
spectrum distortion are expected depending on the type of mass hierarchy (ordering of 
levels) in the scheme: 

a) . The spectrum with the large kink and suppressed tail above the kink. This type of 
distortion realizes in the scheme with inverted mass hierarchy. 

b) . The spectrum with small kink. This type of distortion is expected in the scheme 
with normal mass hierarchy. 

c) . The spectrum with small kink and "strong" bending at the shifted end point. Such 
a distortion corresponds to the non-hierarchical spectrum with normal ordering of states. 

d) . Spectrum with large kink and strongly suppressed tail above the kink and a shift 
of end point. This type of distortion realizes in the non-hierarchical scheme with inverted 
order of states. 

The rates in the suppressed tails and sizes of small kinks are determined by p h e or (1— p£) 
and the latter quantities are restricted by the Bugey or CHOOZ bounds: p h e < 0.027. 
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The lower bound on pi appears from the LSND result. This bound is close to the Bugey 
upper bound in the (3+1) schemes and it can be much weaker for the (2+2) schemes. 

The 4z/-schemes with inverted mass hierarchy are disfavored by the data from SN1987A, 
leading to a hard spectrum of v e . 

5. The identification power of the /3-decay measurements will depend on the possibility 
to detect the suppressed tail and the small-size kinks in the spectrum. 

Even if small kinks or suppressed tails are unobservable, the important conclusions can be 
drawn from comparison of the values of mp and ^ Am 2 LSND . This will allow one to identify 
the type of mass hierarchy (in 4z/-schemes) and also to distinguish the hierarchical from 
non-hierarchical schemes. Note that if the small kink or suppressed tail are unobservable, 
the effect expected from the presently favored schemes of neutrino masses is the same as the 
effect of the electron neutrino with definite mass. 

Observations of the small kinks or suppressed tails will allow us to measure the mixing 
parameters p h e and p[ and to make the independent identification of the scheme. The (5- 
decay measurements can also distinguish the three and four neutrino schemes. However, 
this can be done by MiniBooNE and other experiments even before new /5-decay results will 
be available. 

6. Even if the LSND result is not confirmed, there are some motivations to search for 
the kinks in the energy interval (1 — 10) eV below the end point. The kinks can be due 
to mixing of the active neutrinos with the light singlet fermions which originate from some 
other sectors of theory beyond the Standard model. 

7. The important conclusions can be drawn from the combined analysis of results of 
the /5-decay measurements, 2/30z/-searches and the identification of the solution of the solar 
neutrino problem. 

8. Negative results of future /3-decay experiments will have a number of important 
implications. In particular, 

- Large part of the parameter space of the 3z/-schemes with degenerate mass spectrum 
will be excluded; 



- If the LSND result is confirmed, and the bound from the beta decay is mp < ^Am 2 LSND 

the schemes with inverted hierarchy (order) as well as the schemes with normal order and 



If we want eventually to know the whole "story about neutrinos" we should measure the 
absolute values of their masses. From the point of view of implications for the fundamental 





will be excluded. 
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theory, for astrophysics and cosmology, the masses are at least as important as other neu- 
trino parameters such as the mixing angles and CP-violating phases. As follows from our 
study, to reconstruct the absolute values of masses unambiguously and without additional 
assumptions one needs almost unavoidably to develop and to perform new direct (kinematic) 
measurements (or look for some new alternatives). The KATRIN experiment may be just 
the first step. There is no reason, why we should devote less time, effort and resources for 
determination of the absolute scale of neutrino mass than, e.g., we are going to devote to 
measurements of the oscillation parameters and the CP-violating phases. 



Appendix: Effective mass of the set of quasi-degenerate 
states 

The effect of the set of quasi-degenerate states on the /9-decay spectrum can be described 
by the effective mass, mp, and the coupling p e of the set with the electron neutrino. Let us 
evaluate the accuracy of such an approximation. The error is maximal in the energy interval 
close to Eq — rrii, where rrii is the mass of the lightest state from the quasi-degenerate 
set. Let us compare number of events (decays) due to states from this set in the interval: 
(E — mi — AE) — (E — rrii), using (1) the effective parameters p e and rrip. An(p e , rrip), and 
(2) the precise parameters of states: An(U e i,rrii). Here AE is the energy interval which can 
be resolved by the detector. Let us calculate 

= An(p e , mp) - An(U ej ,mj) 
An(p e ,m /J ) 

Since in the energy range that we are interested, F(E, Z) and p in (0) are smooth functions 
of energy we can factor out them and estimate the relative error as follows 

gj 1^1 2 jggWjo - E)y/(Eo - E) 2 - m 2 Q(E -E- m 3 )dE 

lE °-AE- mi Pe(E - E)^(E -E) 2 -m 2 Q(EQ -E- m p )dE 

Pe [K + AE? - - gj \Uej\ 2 [(m t + AE) 2 - mj]j 
p e [{rrii + AE) 2 - m 2 }l 
where in the sum j runs over the set. 

Note that although the derivative of spectrum at Eq — rrii is divergent, An has finite 
derivative. If \rrij —rrip\ rrip one can expand the relative error around rrip over Am/ ((rrii + 
AE) 2 — ml) as follows 

3 Ej \Ue 3 \ 2 mpAm 3 3 gj \Uej\ 2 [2m 2 - (m r + AE) 2 ](A mj ) 2 
p e [(mi + AE) 2 -m 2 ] + p e [( mi + AE) 2 - m 2 ] 2 ' 1 ' 
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where Arrij = nij — nip. Using the inequality Am <C AE, we can rewrite ([71]) as 

3_Ej l^fm^Am, 3_gj l^fK - 2nijAE - AE 2 ](Anij) 2 

p e (2niiAE + AE 2 ) p e [2 mi AE + AE 2 } 2 ' 1 j 

If bending of the energy spectrum is observable, AE is of the order of m-j or smaller. Con- 
sequently, the first term in ( |7lD is of the order of Am/ AE and the second term is of order 
(Am/ AE) 2 . If we set 







2 




0* 


2 



m/j 

the first term in ([72] ) vanishes and R will be of the order (Am/ AE) 2 . 
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Figure 1: The bounds on the effective /3-decay mass, mp, in the 3z/-scheme with mass 
degeneracy. Shown are the upper bounds from the 2/30^-decay as the functions of mixing 
angle relevant for the solution of the solar neutrino problem (see eq. ( p6|) ). The upper solid 
(dashed) line corresponds to the present bound m ee < 0.34 eV and \U e s\ 2 = (\U e s\ 2 =0.05). 
The lower solid (dashed) line corresponds to m ee < 0.05 eV and \U e 3\ 2 = (\U e ^\ 2 =0.05). 
These lines are drawn in assumption of strong degeneracy of neutrino masses. The vertical 
lines mark the 90 % C.L. borders of the LMA solution region. Shown also are the present 



upper bound on the neutrino mass from structure formation |4l| and the sensitivity limit of 
KATRIN experiment. 
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Figure 2: The four neutrino schemes of mass and mixing. The boxes correspond to the 
mass eigenstates. The position of the box in the vertical axis determines the mass. The 
shadowed parts of the boxes indicate the amount of the electron flavor in the corresponding 
eigenstate that is, |t/ ei | 2 . The solar pair is formed by the two strongly degenerate states 
with Amg and significant amount of the electron flavor. In the hierarchical schemes the 
effective mass of the light set is much smaller than the mass of the heavy set. In the non- 
hierarchical schemes these two masses are comparable. For definiteness we show distribution 
of the electron flavor in the solar pair which corresponds to the large mixing solution of the 
solar neutrino problem. 
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Figure 3: The bounds on mp and p h e in the (3+1) scheme with normal mass hierarchy. The 
dashed curve and solid lines attached to it (from below and above) show the upper bound 
on Pg from Bugey and CHOOZ experiments, respectively. The LSND lower bound (see eq. 
(|3"5j)) is shown by dot-dashed curves. The allowed regions are shadowed. The triplets of 
solid lines show the upper bounds on mp assuming that future 2/30i/-decay searches will give 
m ee < 0.01, 0.03 and 0.05 eV. The central line in each triplet corresponds to the contribution 
from the heaviest mass eigenstate (eq. fl38|)) and the other two lines show the uncertainty 
due to the contribution of light states in the modification of the scheme in which the mass 
hierarchy of the three light states is inverted (see eq. (|39|)). 
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Figure 4: The bounds on rap and p h e in the (2+2) scheme with normal mass hierarchy. The 
dashed curve and solid lines attached to it show the upper bounds on p h e from Bugey and 
CHOOZ experiments, respectively. The LSND lower bound (see eq. (0)) is shown by dash- 
dotted line. The allowed regions of parameters are shadowed. The triplets of solid lines show 
the lower bounds on mp from a positive signal in future 2/50z/-decay searches which would 
correspond to m ee =0.01, 0.03 and 0.05 eV. The central lines correspond to contribution to 
m ee from the heavy set only. 
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Figure 5: The bounds on mp and in the (3+1) scheme with inverted mass hierarchy. The 
dashed curve and solid lines attached to it (from below and above) show the upper bound on 
Pg deduced by Bugey and CHOOZ experiments, respectively. The LSND lower bound (see 
eq. (|3"5j)) is shown by dot-dashed lines. The allowed regions are shadowed. The horizontal 
dashed lines are the upper bounds on m$ from the 2/30^-decay searches which correspond 
to m ee < 0.34 eV, U e3 = 0, and different values of sin 2 20 . 
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Figure 6: The bounds on mp and p h e in the (2+2) scheme with inverted mass hierarchy. 
The dashed curve and solid lines attached to it (from below and above) show the upper 
bound on p h e deduced by Bugey and CHOOZ experiments, respectively. The LSND lower 
bound (see eq. (|35p) is shown by dot-dashed lines. The allowed regions are shadowed. The 
horizontal dashed lines are the upper bounds on mp from the 2/3GV searches which correspond 
to m ee < 0.34 eV and different values of sin 2 26 l . We assume that solar pair of states gives 
the dominant contribution (see eq. fl4T)D). The curve for the SMA solution (not shown here) 
practically coincides with the curve of KATRIN sensitivity. 
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Figure 7: The bounds on the effective mass of the light set mp, and the coupling of the 
electron neutrino with heavy set, p\, in the non- hierarchical schemes with normal order of 
levels. The vertical solid lines show the upper bounds on p^ from Bugey experiment. The 
dashed and dash-dotted vertical lines show lower bounds on p^ from LSND experiment in 
the (3+1) and (2+2) schemes, respectively (see eq. fl35D). The allowed regions for (3+1) 
scheme are shadowed. The lines with different values of sin 2 2# are the upper bounds from 
the 2/?0i/-decay searches which correspond to m ee <0.34 eV. The line denoted by "3+1" 
shows the upper bound for the (3+1) scheme with \U e3 \ 2 =0.05 (see eq. Q5TD), while the 
others are valid both for the (2+2) scheme and the (3+1) scheme with |[/ e 3| 2 =0 (see eqs. 
(|53|), (|51|)). The lines marked by (a) and (b) are calculated for 
0.477 eV, respectively. 
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Figure 8: The bounds on the effective mass of the heavy set, m% and the coupling of 
the electron neutrino with light set, 1 — p%, in the non-hierarchical schemes with inverted 
order of states. The vertical solid lines show the upper bounds on (1 — p^) from Bugey 
experiment. The dashed and dash-dotted vertical lines show lower bounds on (1 — pi) from 
LSND experiment in the (3+1) and (2+2) schemes, respectively (see eq. fl3"5|)). The allowed 
regions for the (3+1) scheme are shadowed. The lines with different values of sin 2 20 are 
the upper bounds from the 2/50z/-decay searches which correspond to m ee <0.34 eV. The 
line denoted by "3+1" corresponds to the (3+1) scheme with |f/ e 3| 2 =0.05 (see eq. |62[) while 
others are valid both for the (2+2) scheme and the (3+1) scheme with |t/ e 3| 2 =0 (see eqs. 

([53]), (|59D). The lines marked by (a) and (b) are calculated for \J Am? LSND = 1.32 eV and 
0.477 eV, respectively. 
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